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Context. Recently, diffuse and extended sources in TeV y-rays as well as in X-rays have been detected in the direction of the Galactic 
globular cluster (GC) Terzan 5. Remarkably, this is among the brightest GCs detected in the GeV regime. The nature of neither the 
TeV nor the diffuse X-ray signal from Terzan 5 is not yet settled. These emissions most likely indicate the presence of several non- 
thermal radiation processes in addition to those that give rise to the GeV signal. 

Aims. The aim of this work is to search for diffuse X-ray emission from all GeV-detected GCs where appropriate X-ray observations 
are available, and to compare the obtained results with the signal detected from Terzan 5. This study will help to determine whether 
Terzan 5 stands out among other GC or whether a whole population of globular clusters feature similar properties. 
Methods. After assessing all archival X-ray observations of GeV detected GCs, we analyzed the data of six Chandra observations 
pointed toward M 62, NGC 6388, NGC 6541, M 28, M 80, and NGC 6139. For each GC we characterized the diffuse X-ray emission 
using the same analysis techniques as for Terzan 5. To study the emission on the same relative scales we used the half-mass radius as 
a scale parameter to determine the extent of the potential emission region. 

Results. None of the six GCs show significant diffuse X-ray emission above the particle and diffuse galactic X-ray background 
components. The derived upper limits allow to assess the validity of different models that were discussed in the interpretation of the 
multi- wavelength data of Terzan 5. A scenario based on synchrotron emission from relativistic leptons provided by the millisecond 
pulsar population cannot be securely rejected if a comparable magnetic field strength as in Terzan 5 is assumed for every GC. However, 
such a scenario seems to be unlikely for NGC 6388, and M62. An inverse-Compton scenario relying on the presence of a putative 
GRB remnant with the same properties as proposed for Terzan 5 can be ruled out for all six GCs. Finally, the assumption that each 
GC hosts a source with the same luminosity as in Terzan 5 is ruled out for all GCs but NGC 6139. 



1. Introduction 

Globular clusters (GCs) are dense systems of very old stellar 
populations. Owing to the extreme densities in their cores, GCs 
feature very high stellar encounter rates that lead to high proba- 
bilities for the dynamical formation of compact binaries ( Pooley" 
|& Hut|[20Q6l IPooley et al. 2003; H ui et al.,,2010 ). This in turn 
likely gives rise to the many of millisecond pulsars (msPSRs) 



Furthermore, a TeV y-ray source was recently detected with 



H.E.S.S. in the vicinity of Terzan 5 (HESS -Collaboration et al. 



that have been discovered in some of these systems ( Camilo & 
|Rasio|2005trRansom|2008| ). 

Recently, several Galactic GCs could be associated to GeV 



y-ray sources detected with the Fermi/LAI observator y ( (Abdo 
[eraL]|2009b| [MTOj |Kong et al.|[2QTQt |Tam et al.|[2QTT] ). The y- 



ray spectra of most of these sources feature cut-offs at a few 
GeV, very reminiscent of the spectra measured from individual 
msPSRs (see, e.g., [Abdo et al.||2009c|a ). Indeed, an individual 
msPSR in a G C has been detect ed recently in the GeV regime 
with Fermi (see Freire et al. 201 1 ). This is one of the reasons why 



the GeV y-ray signals from GCs are generally interpreted as the 
cumulative magnetospheric emission from their whole msPSR 
population. However, the observed emission in the GeV regime 
could also originate from inverse Compton (IC) up- scattering 
of ambient photon fields by leptons accelerated in pulsar wind 



nebulae (PWNe) surrou nding individual msPSRs (Cheng et al 
2010[ Hui et al. 2011| ). In certain cases the resulting spectra 



would look very similar in the GeV range to those observed by 
Fermi from these GCs. 



201 1| ). Intri guingly, this GC is among the brightest GCs seen 
with Fermi (^ Abdo et "STjpOlOl ) and also the one with the most 
detected msPSRs (3^. In the current literature, two models 
predict TeV y-ray emission roughly at a level detected from 
Terzan 5. Similar to what was also suggested for the GeV 
range (see previous paragraph), these models both involve IC 
up- scattering of ambient photon fields, such as starlight from 
GC member stars and the cosmic microwave background, by 
relativistic leptons. As in the GeV regime, these leptons could 
be accelerated either in individual PWNe surrounding msPSRs 
dBednarek & Sitar ek 2007; [Venter et aLl|2009|) or even be re- 
accelerated in colliding PWNe shocks (Cheng et al.|2010 l. 

Leptonic models for TeV y-ray emission from GCs also pre- 
dict X-ray emission caused by synchrotron radiation from the 
same lepton population ( [Venter & de Jagerl|2008| ). Indeed, dif- 
fuse X-ray emission centered on the core and extending well be- 
yond the half-mass-radius (rh) has recently been detected with 
Chandra from Terzan 5 (Eg er et al. 2010 ). The spectrum of this 
emission is well represented by an absorbed powerlaw model 
with a spectral index of 0.9+0.5 and therefore is most likely 
of non-thermal origin. Diffuse X-ray emission has been mea- 
sured from a number of other Galactic GCs (see, e.g., jOkada] 
etaLj2007^. However, in these cases the spectra were softer and 



http://www.naic.edu/~pfreire/GCpsr.html (and references therein) 
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Table 1. The GC sample and Chandra observations 



GC 


Reference*^ 1^ 


h 


ObsID^^) 


Coverage*^"^^ 


Net exposure*^^^ 


Sensitivity*^^^ 


name 




n 




(%) 


(ks) 


(10-15 ergcm-2 s-^) 


M 62 


AlO 


74 


2677 (ACIS-S) 


98 


62.2 


0.6 


NGC 6388 


AlO 


40 


5505 (ACIS-S) 


99 


44.7 


1.0 


NGC 6541 


A10,T11 


71 


3779 (ACIS-S) 


80 


44.9 


0.7 


M 28 


AlO 


94 


9132 (ACIS-S) 


63 


141.9 


0.5 


M 80 


Til 


37 


1007 (ACIS-S) 


99 


48.6 


1.0 


NGC 6139 


Til 


51 


8965 (ACIS-S) 


83 


18.1 


2.0 


47Tuc 


A09,A10 


190 


Th too large 








0) Cen 


AlO 


134 


Th too large 








NGC 6440 


AlO 


35 


bright point source 








NGC 6441 


A10,T11 


34 


bright point source 








NGC 6624 


Til 


49 


only grating or low exposure 








NGC 6652 


AlO 


29 


only low exposure 








M 15 


AlO 


60 


bright point source 








NGC 6752 


A10,T11 


115 


Th large, known diffuse source in region ( Okada et al. 2007 1 




Liller 1 


Til 


48 


only grating 









Summary of GeV-bright GCs that were conside red (t op section) or not considered (bottom section) in this work. ^^"^ Fermi/LAI reference: A09 
( |Abdo et al.|20 09b ), AlO jAbdo etal.|20Tot , Til ( |Tam et al. 201 1 ) (2)half-mass-radius (rh) from |Harris|fl996, 2010 edition| ). Chandra observation 
identifier and instrument mode, or reason for exclusion. *^^^Fraction of the area between one and three within t he Fo V. '^^^emaining exposure after 
GTI screening. Owing to low background activity this is identical to the total exposure in most cases (see Sect. |2.2| ). Sensitivity for point-source 
detection (0.5-7.0 keV). 



the emission region was either offset from the core in the direc- 
tion of the GCs proper motion or featured an arc-like morphol- 
ogy. Here, the emission was interpreted as either non-thermal 
bremsstrahlung from shock-accelerated electrons hitting nearby 
gas clouds, or thermal e mission from shoc k-heated gas ( Okada 
et aI.||2QQ7] ). Following Eger et al.| ( |2Q1Q| ), these scenarios are 
most likely not valid for Terzan 5, whereas synchrotron radiation 
appears to be a plausible explanation considering the energetics. 
Based on a study of available multi- wavelength data of Terzan 5, 
predominantly in the radio regime, [Clapson et al.| ( |2Q1 1 ) were 
able to disfavor non-thermal bremsstrahlung, while still leaving 
other non-thermal processes, such as IC or synchrotron radia- 
tion, as valid scenarios. Following " Cheng et al.| ( |2010 ), the dif- 
fuse X-ray signal could be interpreted as the tail of the IC peak 
at TeV energies. However, these authors' prediction for the flux 
is a factor of 3-4 lower than what is observed from Terzan 5. 

Another possible explanation for the TeV gamma-ray emis- 
sion is related to a remnant of a short gamma-ray burst (ORB) 
(Domainko 201 1)) resul ting from a neutron star-neutron star 
merger ( Grindlay et al . |[2006] ) . Owing to their high percentage 
of stellar binaries, GCs are prime candidates to host these events 
because they form considerable numbers of potential progeni- 
tors. Here the y-rays could be produced in inelastic collisions 
of hadrons, accelerated in the shocks, with ambient interstellar 
matter. Then the TeV y-ray signal could be accompani ed by ther- 
mal X-ray emission arising from shock-heated gas (Domainko 
i& RuflTert 2 008) or non-thermal IC radiation from primary lep- 
tons accelerated along with the hadrons in the shocks ( Domainko] 
201 1| ). Particularly for the latter case the TeV and X-ray signals 
that have been detected from Terzan 5 could be related. 

The aim of this work is to systematically search for diffuse 
hard X-ray emission from other Galactic GCs that have been 
detected in the GeV energy range by Fermi. Based on differ- 
ent model assumptions we then investigate whether this kind of 
emission could be typical for a whole population of GCs or if 
Terzan 5 may stand out, as this GC would do if there were a 
recent and rare catastrophic event. To facilitate a comparison of 



the obtained results to Terzan 5, we studied the diffuse emission 
on similar relative scales, using rh as a reference parameter to 
determine the size of each potential emission region. 



2. X-ray observations, analysis and results 

2.1. The globular cluster sample and Chandra observations 

We built the list of GeV-bright GCs d etected by Fermi b ased 
on the systematic studies performed by Abdo et al. ( 2010 ) and 
Tam et alT] ( |2011| ), regardless of the measured spectral proper- 
ties and their interpretation. Furthermore, we restricted our study 
to Chandra data ( Weisskopf et al. 2002) because it is the only 
X-ray telescope that can cope with the task of effectively re- 
moving point-like X-ray sources in crowded areas close to GCs. 
The narrow point spread function (PSF) of Chandra's telescope 
keeps the potential contamination of regions outside the removed 
point sources to a minimum. Fortunately, all selected GCs have 
been observed by Chandra and the data were accessible via the 
public archive. To retain as much spatial and spectral informa- 
tion as possible at the same time, we only considered the non- 
grating Adv anced CCD Imaging Spectrometer (ACIS, jGarmire 
et al.|2003] ) data with an exposure of at least 10 ks. 

To judge which of these objects are suited for our study we 
required that more than 50% of the region of a potential dif- 
fuse X-ray signal is covered by the ACIS field of view (FoV). 
For Terzan 5 the significant diffuse excess emission lies between 
one and three rh, a scale we used to define the coverage frac- 
tion for all other GCs. The respective values for r^ were taken 
from Harris ( 1996, 2010 edition). Furthermore, we required that 
no X -ray point source brighter than 1x10"^^ erg cm"^ s"^ (see 
Sect. |2.2| ) is located within the FoV to avoid contamination from 
the PSF wings in our region of interest. For NGC 6752 and M 80, 
already known dif fuse X-ray sources are located within 3xrh 
( Okada et al.|2QQ7) , which would have to be removed from the 
data to search for an additional signal. However, for NGC 6752 
the area of the remaining region is below the required minimum 
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coverage fraction of 50% and this GC was therefore excluded 
from the analysis. The sample of GeV-bright GCs is compiled in 
Table [T] where the top section lists all GCs passing our a priori 
cuts and the bottom section gives short comments on why the 
respective GC was excluded from our study. 

2.2. Data preparation and point-source detection 

All steps of the analysis procedure were performed in the same 
way as for our study of Terzan 5 (see 'Eg er et aT]|2010| . For 
the X-ray analysis we used the CIAO software version 4.3, sup- 
ported by tools from the FTOOLS package (Blackburn 1995) 
and XSPEC version 12.5.0 for spectral modeling ( Amaud,1996 ). 
The event 1 data were reprocessed with the latest position- and 
energy calibration (CTI correction, v4.4.6) using bad pixel files 
generated by acis_runJiotpix. To screen the ACIS-S data 
from periods of increased background flaring activity, we used 
the good-time-interval (GTI) file that is provided by the standard 
processing chain. In most of the observations the background 
level was fairly low, and the net exposures after the screening 
are virtually identical to the total lifetime. To test the standard 
GTI filtering we inspected the light curves created from events 
in the 9-12 keV band for each observations. However, no ad- 
ditions to the standard GTIs were necessary for any of the six 
datasets. The resulting net exposures are also listed in Table [T 

To detect and remove point-like X-ray sources, we ran the 
wavelet detection algorithm wavdetect in the three energy 
bands 0.5-2 keV, 2-7 keV, and 0.5-7 keV. These bands were 
chosen to be sensitive to soft Galactic sources, such as active 
stars, hard sources associated to the GCs, and absorbed back- 
ground sources, e.g. AGN. The sensitivities for the detection 
of point sources are also given in Table [T] To estimate whether 
the detected sources are predominantly related to the GC or be- 
long to a population of background sources, such as AGN, we 
used the results from [Giacconi et al.| ( |2001| ) to calculate the 
estimated number density of background sources based on the 
detection sensitivity of the individual observations. We found 
that the number densities of expected background sources lie 
between -1900 sources deg"^ (M 62) and 675 sources deg"^ 
(NGC 6139). For an individual ACIS CCD chip this yields total 
numbers of expected sources of 34 (M 62) and 12 (NGC 6139), 
respectively. These results are compatible with similar calcula- 
tions for Terzan 5 made by Heink e et al.| ( |2006| , who expected 
a total number of eight AGN per CCD chip above a sensitiv- 
ity threshold of ~3xlO"^^ erg cm"^ s"^ (^^10 counts). We note 
that the galactic column densities in the directions of the GCs 
considered here are about an order of magnitude lower than to- 
ward Terzan 5. Therefore, we did not correct the number of ex- 
pected background sources for absorption eff'ects. Except for the 
region within rh these estimated number densities of background 
sources agree very well with the number of detected sources, im- 
plying that most of these sources are of extra-galactic origin. 

To estimate the fluxes of these sources based on their count- 
rate, we assumed a typical non-the rmal spectrum, suc h as mea- 
sured for X-ray sources in Terzan 5 ( [Heinke et al.|2006' ). Here we 
assumed an intrinsic powerlaw spectrum with index -2 and took 
into acc ount the total Galactic column density (A^h, from Dickey 
& Lockman 1990) in the respective direction. Among all con- 
sidered GCs the values of A^h are moderate and on the order of a 
few times 10^^ cm"^. To estimate the flux, we first corrected the 
count-rates for the eff'ects of the mirror vignetting based on the 
off'-axis angle at each source position. If the derived flux of any 
of the point sources in the FoV was above a level of 1x10"^^ erg 
cm"^ s"^ the respective GC was excluded from the sample (see 



Sect. |2.1| ). For the remaining six GCs we removed all detected 
point sources from the data using the 3 cr width of the PSF at the 
source position. 

As already mentioned in Sect. |2.1[ an al ready known diff'us e 
X-ray source is located close to M80 (see Qkada et"aL]|2007 ). 
This source is located toward the southwest of the GC center 
and the significant emission extends out to ~4rh. To remove this 
source from the data we excluded a circle segment with an open- 
ing angle of 90° oriented toward the southwest (270°-360°, west 
to north) and extending out to 5rh. 

2.3. Diffuse emission 



To characterize the level and radial dependence of potential dif- 
fuse X-ray emission from each of the six GCs on the same rel- 
ative scales as for Terzan 5, we extracted spectra from the same 
eight annular regions ("Ringl" to "RingS") with inner and outer 
radii in units of rh. Because the X-ray point source populations 
of the GCs are much smaller compared to Terzan 5, we also con- 
sidered regions closer than Ixrh that had no high contamination 
from unresolved point sources. Therefore, we introduced one ad- 
ditional annulus ("RingO") closer toward the GC center with the 
same widths as the other eight annuli. We chose a scaling of 
the region sizes with respect to r^ because it reflects the typi- 
cal size of the potential emission region for all scenarios. For IC 
emission from relativistic leptons the emission is expected to fol- 
low the profile of the target photon field, whereas in an msPSRs 
scenario the emission should reflect the distribution of sources 
within the GC. We used the values for the position and rh as 
given in the GC catalog of Harris|( |1996, 2010 editionj ). Table |2] 
lists the details of these regions. The mean eff'ective area (ARE) 
and energy resolution (RMF) files were calculated by averaging 
over all contributing pixels assuming a flat detector map. The 
AREs and RMFs created using a weighted average over all pixels 
with weights corresponding to acceptance-corrected count-rates 
yielded comparable results. 

While the position of the GC core is always located in the 
area covered by the back-illuminated ACIS-S3 chip, in some 
cases parts or all of outer three annuli were on one of the two 
neighboring front-illuminated chips. In these cases we extracted 
the spectrum and response files for each CCD separately to fit 
the datasets in parallel (see next paragraph). As for Terzan 5 (see 
Eger et al. 2010 ), to estimate the non-X-ray background (NXB) 
component we used a background dataset from the calibration 
database where the ACIS detector was operated in its stowed po- 
sition. For each spectrum from the annular regions a correspond- 
ing background spectrum was extracted from the NXB dataset 
also using the identical point- source exclusion regions. To ac- 
count for the time-dependence of the NXB, we scaled the back- 
ground by the ratio of the source and background count-rates in 
the 9-12 keV energy band for each spectrum (as described by 
Markevitch et al. 2003). 

We fitted the NXB -subtracted spectrum from each annulus 
with an absorbed powerlaw model. Owing to lack of statistics in 
individual regions, the goal was only to reproduce the observed 
spectra reasonably well with a simple model to derive observed 
fluxes. For each GC we linked the hydrogen column density A^h 
for all eight rings, but left the photon index and the normaliza- 
tion free to vary independently. If more than one spectrum was 
available for a given annulus due to two or three contributing de- 
tector CCDs, we linked all spectral parameters, with a fixed ratio 
between the normalizations, corresponding to the relative size of 
the extraction areas on each CCD. As for Terzan 5, because of 
the limited statistics in the spectra from individual annuli, the ab- 
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Table 2. Parameters of the annular extraction regions 



Table 3. Spectral fitting results 





M62 


NGC 


NGC 


M28 


M80 


NGC 






6388 


6541 






6139 


width^i^ 


36 


19 


34 


45 


18 


25 


radiusO^^^ 


46 


25 


44 


58 


23 


32 


radius 1 


81 


44 


78 


103 


40 


56 


radius2 


117 


63 


112 


148 


58 


81 


radiusS 


152 


82 


146 


193 


75 


105 


radius4 


186 


101 


179 


236 


92 


129 


radius5 


222 


120 


213 


281 


110 


153 


radius6 


258 


139 


247 


326 


127 


177 


radius? 


293 


158 


281 


371 


145 


202 


radius 8 


329 


178 


315 


416 


163 


226 



All values are given in arc seconds. ^^^Width of the annuli equals 
0.48xrh for each GC. *^^^Inner radius of each ring calculated as (0.62, 
1.10, 1.58, 2.06, 2.52, 3.00, 3.48, 4.00, 4.44) xr^. 



sorbed powerlaw model reproduced the spectra suflftciently well 
with reduced values in the range of 1.1 to 1.4. We calculated 
the surface flux for each region by dividing the measured total 
flux by the efl'ective extraction area, taking into account excluded 
regions, regions outside the FoV as well as bad columns and pix- 
els. Table[3]contains the spectral fitting results for those four GCs 
where significant excess emission above the NXB was detected. 
The radial evolution of the surface brightness is shown in Fig.[T] 
Here we show the observed fluxes, because the absorbed power- 
law model is only a rough approximation of the spectral shape. 
A correction for absorption would require a fit of individual an- 
nuli with a more physically motivated model, which the limited 
statistics does not allow for. However, assuming that the fore- 
ground column density does not vary much within the regions, 
the observed and intrinsic fluxes can be expected to scale simi- 
larly. 

The radial difl'use surface flux of M 80 seems to be peaked 
in the fourth bin at about IOC. We investigated this by visual 
examination of the raw images and smoothed versions in difl'er- 
ent energy bands for any point-like or extended but still compact 
source located in this annulus that might have escaped the detec- 
tion algorithm somehow, but we did not find any. However, de- 
spite this diff'erence, the fluxes in all annuli are still comparable 
within the uncertainties. Possible explanations for this increase 
are statistical fluctuations or a faint extended source, such as a 
background cluster of galaxies with a flux below the detection 
threshold. 

To investigate the nature of the observed difl'use X-ray emis- 
sion in more detail, we fitted the combined spectrum extracted 
from the outer three regions with a more physically motivated 
model. In Terzan 5 we found significant diff'use X-ray excess 
emission above the Galactic diff'use background only within the 
inner five annuli, whereas the combined spectrum from the outer 
three annuli was compatible with pure diff'use Galactic back- 
ground. Even though we did not detect any significant depen- 
dence of the surface flux on the distance to the GC center, we 
still considered the inner five annuli as our source region of po- 
tential excess emission and the outer three annuli as the region 
containing predominantly background emission. 

Here we performed the same study as for Terzan 5 in us- 
ing a non-equilibrium ionization model (NEI in XSPEC) to de- 
scribe the emission from diff'use hot Galactic gas. This model 
was already used by [Ebisawa et al.| p005[ henceforth EOS), 
who performed deep Chandra observations of a typical "empty" 



Region 


Excess counts^^^ 




Surface flux*^^^ 


fd f ) 


M62 










RingO 


660.1 +43.3 


8.9 + 0.3 


2.51+0.42 


1.2(504) 


Ringl 


897.4 + 49.5 


9.6 + 0.1 


2.18 + 0.26 




Ring2 


1145.9 + 48.1 


9.8 + 0.1 


2.05 + 0.20 




Ring3 


1175.2 + 50.8 


9.9 + 0.1 


1.82 + 0.12 




Ring4 


1669.7 + 58.3 


9.6 + 0.2 


2.22 + 0.19 




Ring5 


1750.0 + 78.7 


9.7 + 0.1 


2.18 + 0.18 




Ring6 


850.1 +41.8 


9.0 + 0.4 


1.91 + 0.15 




Ring7 


223.9 + 14.8 


8.6 + 0.5 


2.18 + 0.25 




Ring8 


326.5 + 17.7 


9.0 + 0.6 


2.13 + 0.20 




Outer 


1672.5 + 74.8 


0.40 + 0.03 


5.74 + 0.37 


1.3(99) 


NGC 6388 








RingO 


264.9 + 26.5 


5.6 + 1.2 


6.65 + 2.4 


1.1(469) 


Ringl 


318.0 + 37.8 


6.8 + 1.1 


3.79 + 1.3 




Ring2 


392.2 + 38.1 


7.0 + 0.9 


3.34 + 0.91 




Ring3 


446.8 + 37.0 


6.6 + 0.9 


3.72 + 0.98 




Ring4 


457.6 + 29.6 


7.1 + 0.9 


2.63 + 0.66 




Ring5 


580.4 + 46.2 


6.9 + 0.9 


3.03 + 0.71 




Ring6 


681.4 + 46.6 


7.4 + 0.9 


2.76 + 0.61 




Ring7 


705.5 + 55.4 


7.1 + 0.8 


3.20 + 0.64 




Ring8 


657.7 + 45.6 


7.3 + 0.9 


O CO 1 f\ cc\ 

2.58 + 0.59 




Outer 


1966.9 + 76.5 


0.18 + 0.02 


8.25 + 0.55 


1.1(88) 


NGC 6541 








RingO 


442.5 + 47.3 


8.9 + 0.3 


1.56 + 0.19 


1.4(451) 


Ringl 


730.2 + 53.8 


9.2 + 0.2 


1.42 + 0.13 




Ring2 


917.3 + 63.4 


9.5 + 0.1 


1.34 + 0.10 




Ring3 


769.7 + 58.8 


9.2 + 0.2 


1.39 + 0.12 




Ring4 


437.4 + 21.0 


9.5 + 0.1 


1 or\ 1 f\ f\c\ 

1.20 + 0.09 




Ring5 


547.5 + 23.1 


9.1 + 0.1 


\A1 ± 0.10 




Ring6 


450.1 +21.6 


9.2 + 0.2 


1.38 + 0.11 




Ring7 


367.1 +20.0 


9.4 + 0.2 


1.25 + 0.11 




Ring8 


345.0+ 18.9 


9.4 + 0.2 


1.25 + 0.11 




Outer 


1767.8 + 81.0 


0.17 + 0.02 


5.51 + 0.39 


0.9(67) 


M80 










RingO 


187.6 + 24.6 


2.5 + 1.6 


10.8 + 8.0 


1.4(567) 


Ringl 


210.6 + 28.5 


2.6 + 1.4 


6.9 + 5.7 




Ring2 


291.6 + 28.4 


3.3 + 0.6 


3.1 + 1.5 




Ring3 


360.6 + 32.2 


3.2 + 0.5 


3.9 + 1.7 




Ring4 


432.2 + 35.1 


1.6 + 0.9 


16 + 11 




Ring5 


538.8 + 42.5 


2.2 + 0.5 


8.1+3.6 




Ring6 


626.3 + 43.2 


2.0 + 0.6 


10 + 5.0 




Ring7 


673.4 + 51.3 


2.5 + 0.4 


6.0 + 2.2 




Ring8 


711.8 + 45.5 


2.2 + 0.4 


9.4 + 3.4 




Outer 


2309.4 + 79.0 


0.42 + 0.08 


26.0 + 2.5 


1.2(182) 



The quoted uncertainties are 90% confidence intervals calculated by 
letting each free parameter vary. Therefore, parameter correlations are 
taken into account. *^^^Total number of excess counts summed over all 
contributing CCDs in the energy bands 1-7 keV (Ring 1-8) and 0.7-10.0 
keV (Outer). ^^^This column gives the photon index in the cases where 
an absorbed powerlaw spectrum was fitted (Ring 1-8), or kT in units of 
keV for the combined Outer region. Observed surface flux in units of 
10~^ ergcm"^ s~^ sr~^ in the energy bands 1-7 keV (Ringl-8) and 0.7- 
10.0 keV (Outer). 



Galactic Plane region and in particular studied the diff'use emis- 
sion after the removal of point sources in great detail. In con- 
trast to the powerlaw fitting of individual annuli here we chose 
an energy interval of 0.5-10 keV to be able to directly compare 
our results to those obtained by EOS. The column density and 
consequently also the expected amount of emitting gas is about 
an order of magnitude higher toward the region studied by EOS 
compared to the observations used here (A^h = 2.07x10^^ cm"^ 
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vs. A^H ~ 1-2x10^^ cm"^, respectively). To account for the lack 
of statistics arising from shorter exposures and lower fluxes, we 
fixed most of the model parameters to the values given in Table 8 
from EOS by leaving only the normalization and the temperature 
free to vary. To fit the spectra from the outer regions, we first 
started with only one temperature component and that already 
yielded good results for three GCs (M 62, NGC 6388, NGC 
6541). In these cases the derived temperatures are comparable 
to the soft component in EOS. For M 80 a second component 
was needed to describe the data and we used the hard component 
from EOS with a temperature fixed at their value {kT = 5 keV). 
To compare the integrated observed surface fluxes to EOS, we 
chose the same energy range (0.7-10.0 keV). These results are 
also listed in in Table [S] As discussed in detail in sect. |3.H we 
conclude that the difl'use X-ray emission in all cases is compat- 
ible with the difl'use Galactic emission and no extra component 
above this background is found. 



Table 4. GC upper limits 



GC 


Counts '^^^ 


Excess Rate'^^ 


T TT Z7(3) 

UL F\,' 


name 




(lU s ) 


(10 erg cm s ) 


Rings 1-5 








M62 


13392 


9.29 


1.1 


NGC 6388 


3749 


6.85 


0.89 


NGC 6541 


6681 


9.10 


1.1 


M28 


49219 


7.80 


0.96 


M80 


2591 


5.24 


0.63 


NGC 6139 


2197 


13.0 


1.7 


Rings 0-5 








M62 


14502 


9.67 


1.2 


NGC 6388 


4121 


7.18 


0.93 


NGC 6541 


7385 


9.60 


1.2 


M28 


55262 


8.28 


1.0 


M80 


2876 


5.52 


0.66 


NGC 6139 


2386 


13.5 


1.7 



2.4. Upper limits 



None of the six GCs shows signs of significant diff'use excess 
emission between 0.6 and 4.4 above the NXB and Galactic 
diff'use background components. As mentioned in Sect. [T] this 
region could be interesting, because a significant diff'use signal 
was detected from the GC Terzan S ( |Eger et al.^2 010). To derive 
upper limits for diff'use excess emission from each GC, we used 
the total number of counts in the 1-7 keV band within the re- 
gions Ringl-S and RingO-S, respectively, again using the same 
exclusion regions as for the spectra. The first region was cho- 
sen such that the derived upper limits can be directly compared 
to models based on the level of diff'use X-ray excess flux from 
Terzan S, which was measured from the respective region. The 
latter area (RingO-S) also takes into account the additional re- 
gion of RingO, which was not possible for Terzan S because of 
its large point source population compared to the GCs studied 
here. The upper limit from the RingO-S region would allow one 
to constrain models for difl'use X-ray flux that are not based on 
the Terzan S measurement, but only rely on other already known 
parameters of the respective GC. Because in many cases the dif- 
fuse X-ray flux is expected to drop rapidly with increasing dis- 
tance from the GC center, the additional area of RingO, which is 
located closer to the GC cores, would produce tighter constraints 
on these models. 

Assuming that all the counts in the respective regions are 
due to background, we calculated the required number of ex- 
cess counts to yield a detection with a statistical significance of 

5 cr. We converted the excess count-rate into a flux by assum- 
ing an absorbed powerlaw model with an A^h value from [Dickey | 

6 Lockman] ( |1990| for each GC position and a photon index of 
-0.9 (as measured for Terzan S, see E ger et al.||2010| ). To cal- 
culate flux upper limits for the total geometric area of RingsO-S 
and Rings 1-S we scaled the values with the inverse of the cover- 
age fraction, assuming a homogeneous flux distribution over the 
whole region. Note that the coverage fractions quoted in Table [T] 
also include excluded regions from point-like or difl'use sources 
as well as bad pixels and columns. Furthermore, we converted 
the flux upper limits into upper limits on the luminosity using 



Total number of counts in the inner five (top section) and inner six 
(bottom section) rings, respectively. *^^^Required excess rate for a 5 cr 
detection. ^^^Flux upper limit (1-7 keV). Note that these upper limits on 
the flux are not yet corrected for the coverage fraction, but the luminos- 
ity upper limits quoted in Table [5] are. 



3. Discussion 

In this section the implications of the above analysis of difl'use 
X-ray emission from the six GCs is discussed. First, we com- 
pare the observed surface fluxes and spectra from the four GC 
where such an analysis was possible to what might be expected 
from pure Galactic difl'use emission from hot gas. The rest of 
the discussion will then evaluate three models for difl'use X- 
ray excess emission related to the GCs in the context of the 
derived upper limits. We also note here that the population of 
unresolved point sources below the detection threshold of the 
Chandra observations would contribute to any extended X-ray 
emission. Therefore, the quoted upper limits for difl'use emission 
may be considered as quite conservative in this sense. Because 
we particularly aim to determine whether the X-ray signal de- 
tected from Terzan S is a general feature of GCs or whether 
it stands out as an exception, we use the flux measured from 
Terzan S as a reference. Because there is no measurement of 
difl'use X-ray flux from Terzan S corresponding to RingO, we 
only used the upper limits derived from the region covered by 
Rings 1-S to compare with the model predictions in sections 3.2 
toK4\ 



the distances listed in the GC catalog of Harris ( 1996, 2010 edi- 
|tion| ). Table|4]contains the information on count statistics and the 
derived upper limits for each of the six GCs. 



3.1. Galactic diffuse emission 

Before the discussion of potential excess emission above the var- 
ious background components, we first examine the NXB and 
diff'use Galactic background. Using the GTI filtering and the 
"stowed" dataset, we are confident that the NXB has been re- 
moved as best as possible from the spectra. Instead of removing 
the Galactic diff'use background as well, we decided to model the 
emission because this component is much more dependent on 
the actual observation position and might vary in spectral shape 
as well as in intensity between diff'erent fields. 

In principle, one could compare the fluxes and tempera- 
tures derived from the NEI model fit with dedicated studies of 
Galactic diff'use emission, such as performed with Chandra by 
EOS. As already mentioned in Sect.|2.2| the column densities to- 
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CO 
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core distance (arcsec) 
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50 100 

core distance (arcsec) 
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Fig. 1. Radial dependence of the surface flux (1-7 keV) above the NXB for four of the six GCs considered in this study (as labeled 
in the individual figures). The diff'erent colors and markers show the measured surface flux (crosses, black) and the predictions 
based on an IC model (squares, magenta), a model assuming the same luminosity as Terzan 5 (stars, red) and a model assuming a 
connection to msPSRs (triangles, blue). The drawn uncertainties are the same as given in Table [3] 



ward the regions here and the field studied by EOS are about a 
factor of ten lower. Assuming that the amount of X-ray emitting 
gas is proportional to the amount of total interstellar material, 
we scaled the flux measured by EOS with the ratio of A^h- For all 
but one GCs (M 80) the derived expected flux is about a factor 
of ~S below the measured value. This could be explained by the 
fact that most of the matter is also contributing to the absorption 
of the X-ray flux, which is strongly dependent on the spectral 
shape of the intrinsic spectrum and the line-of-sight geometry. 
Particularly for soft spectra (kT < 1 keV) absorption has a very 
strong impact on the observed flux and it scales highly non-linear 
with A^H- We tested the relative influence of the intrinsic flux 
and the foreground absorption on the observed flux for our case 
based on the model parameters from EOS. Increasing both A^h 
(from 10^^ cm"^ to 10^^ cm"^) and the intrinsic flux by a factor 
of ten the observed flux from the soft component (kT = 03 keV) 
even decreases by 26%, whereas the observed flux from the hard 
component (kT = 5 keV) scales nearly as its intrinsic flux. This 
means that the flux expectation based on a scaling with Nu is also 
highly dependent on the relative contribution from both compo- 
nents. 

We found that for three of the four GCs (M 62, NGC 6388, 
NGC 6S41) the expected surface flux is compatible with the re- 



sults from EOS taking the efl'ects of absorption into account. 
Because a single (soft) component was sufficient to reproduce 
the spectra from M62 and NGC 6388, we also only took this 
component into account. For M 80, however, the observed sur- 
face flux is at the same level as found by EOS for a region with 
17 times the A^h- This could be explained by the hard component 
contributing more than in the region studied by EOS, which is 
also supported by the fact that only for this GC the addition of a 
S keV thermal component greatly increased the quality of the fit. 
However, the data for M 80 have a much lower statistical quality 
compared to the region studied by EOS because of the shorter ex- 
posure time. Consequently, a more detailed study of the relative 
contribution and separate foreground absorption of the two com- 
ponents is not possible with the available data. Als o, the diff'use 
thermal source in the vicinity of M 80 detected by Qkada et al. 
( [2007| ) might be more extended than suggested by their image 



and could therefore contaminate our region of interest. In this 
case as well, only a dataset with a significantly longer exposure 
would be needed to better determine the extent of this relatively 
weak diffuse source. 
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Table 5. GC model predictions 



3.3. Scenario 2: IC emission from a GRB remnant 





^(1) 
a 




TTT 7^(3) 


Tr(4 


mSr oKS^ 


name 


(kpc) 




(10^3 


erg s~^) 




M62 


6.8 


0.62 


0.57 


15+2 


0.81+^-5^ 

-0.38 


NGC 6388 


9.9 


1.1 


0.99 


19+3 


1 9+2-^ 


NGC 6541 


7.5 


0.99 


0.89 


8.3+1.2 


27+0-2^ 


M28 


5.5 


0.57 


0.51 


6.0+0.9 


46+0'^^ 


M80 


10.0 


0.67 


0.61 


6.4+0.9 




NGC 6139 


10.1 


2.3 


2.1 


7.2+1.0 





( 11996, 2010 ed ition). ^^^Luminosity 
^^^Luminosity upper limit (1-7 keV) 



*^^^Distance to the GC from Harris \ 

upper limit (1-7 keV) for rings 0-5. "^^^Luminosity upper limit (1 
for rings 1-5. "^"^^ Luminosity prediction for an inverse Compton scenario. 
Luminosity prediction for a scenario based on millisecond pulsars. 



3.2. Scenario 1: Synchrotron radiation from leptons 
produced by msPSRs 

The TeV y-ray signal detected from the direction of Terzan 5 
may be explained by IC emission from leptons accelerated either 
in pulsar wind nebulae (PWNe) surrounding individual msPSRs 
([Bednarek & Sitarek|2QQ7|) or in colliding PWN s hocks ( |Cheng| 
et al.||2Q10| ). Following [Venter & de Jager| ( [2QQ8 l), the IC emis- 
sion in the TeV band is supposed to be accompanied by X-ray 
emission arising from synchrotron radiation from the same lep- 
ton population. However, in simple models the apparent offset 
of the TeV signal from the core of Terzan 5 challenges this inter- 
pretation, because the IC signal would be expected to follow the 
morphology of the GC radiation field as well as that of the lepton 
sources, which would also be traced by the X-ray emission. 

Nevertheless, assuming leptons generated by msPSRs are 
indeed responsible for the X-ray and TeV emission seen from 
Terzan 5, this scenario can be tested for all six GCs. If the mag- 
netic field strength is similar among all GCs, the X-ray luminos- 
ity then only scales with the number of msPSRs. The total num- 
ber of msPSRs can be derived by assuming that the GeV signal 
detected by Fermi/ LAT is caused by magneto spheric emission 
from individual pulsars ( |Abdo et al.|2Q10| ). The total number of 
msPSRs then is proportional to the luminosity in the GeV range. 
Therefore, to estimate the X-ray luminosity for each of the six 
GCs we scaled the luminosity measured from Terzan 5 with the 
ratio of the GeV luminosities. For this we used the results ob- 
tained by|Abdo et al.| ( |20TQ| ) (M62, NGC 6388, M28) and |TSn| 
leFaTlpOllI ) (NGC 6541, M 80, NGC 6139), respectively. Our 
luminosity predictions can be found in Table [5] (column labeled 
"msPSRs"). The confidence intervals given in the table are based 
on the uncertainties of the Fermi and Chandra measurements. 

Only for the two GCs M 62 and NGC 6388 the model predic- 
tions lie clearly above the upper limits (UL15). However, because 
of the relatively large error intervals of the Fermi measurements 
the values barely agree within the uncertainties. Therefore, we 
conclude that we cannot rule out a model based on msPSRs for 
any of our GCs, even though there are indications for discrepan- 
cies between upper limit and prediction for two of them. Because 
the dominant uncertainties come from the Fermi measurements, 
improved constraints on a msPSR scenario can be achieved from 
a GeV dataset with better statistics. 



The TeV y-ray signal from Terzan 5 might be associated to a 
recent short GRB event resulting from a neutron stars merger 
(as discussed in Domainko 2011). In this scenario the TeV y- 
rays would be produced in inelastic collisions of hadronic cos- 
mic rays accelerated in the relativistic shock waves with ambient 
interstellar matter. Here the hard diffuse X-ray emission could 
arise from IC up- scattering by already cooled and mildly rela- 
tivistic primary leptons accelerated at the GRB shock waves, if 
the lepton to hadron ratio is about 0. 1 . 

Assuming that such a short GRB remnant is also present in 
our GCs, we calculated the expected X-ray flux based on the 
same lepton energy density required to explain the emission 
from Terzan 5. In this case, the X-ray luminosity would only 
scale with the energy density of the target photon field, where 
for GCs the main contribution comes from the starlight from 
their densely populated core regions. To calculate the expected 
X-ray luminosities for each of the six GCs, we converted the to- 
tal measured diffuse X-ray flux of (5.5+0.8)xlO~^^ erg cm"^ s"^ 



(1-7 keV) from Terzan 5 ( Eger et al.|2010| ) into a l uminosity us- 
ing the updated distance estimate of 5.9 kpc from ( [Ferraro et al. 
[2009) and scaled it with the rat io of the optical luminosity given 
by [Harris] ( [1996, 2010 edition ). The resulting X-ray luminosity 
for this GRB/IC scenario are given in Table [5] (column labeled 
"IC"). The uncertainties here are based on the 1 cr statistical un- 
certainty of the flux measurement from Terzan 5 with Chandra. 
We also show the expected radial surface flux profile based on 
the measurements from Terzan 5. For this purpose we performed 
the same scaling as for the total source on the luminosities con- 
tained in each individual annulus. Then we calculated the ex- 
pected surface flux from each annulus taking the distance to each 
GC and the surface area of the ring into account. For this model 
the predictions are shown as squares (magenta) data points in 
Fig.[T] Errors only reflect the uncertainty of the measured radial 
profile from Terzan 5. 

It becomes clear from the values that an IC emission scenario 
from a putative GRB remnant with the above assumptions can be 
ruled out for every GC. The high values for the expected X-ray 
luminosities arise from the fact that all GC exhibit higher optical 
luminosities than Terzan 5. That leads to more efficient IC up- 
scattering in these systems. 

Two conclusions can be drawn here. Firstly, the non- 
detection of diffuse X-ray emission in our GC sample could 
result from the low rate of GRBs in GCs. As has been dis- 
cussed in Domainko, {201 1 ^ and references therein), the rate of 
short GRBs in the Milky Way is (0)lO-^(f-^ /lOO)-^ per year, 
where 1 ^ fh < 100 is the beaming factor of short bursts. 
Consequently, the probability to find two young GRB remnants 
in the Milky Way is low. Secondly, if extended TeV emission in 
another GC were to be detected in the future and this emission 
can again be interpreted as a GRB remnant, this would point 
toward different remnant properties in the two systems. Gamma- 
ray burst remnants could differ in the total energetics or the pro- 
ton to electron ratio. 



3.4. Scenario 3: A source with the same iuminosity as 
Terzan 5 

The perhaps simplest assumption is that a source with the same 
X-ray luminosity as Terzan 5 resides in each GC. Furthermore, 
in contrast to the previous scenario we assume that the X-ray 
signal does not arise from IC up- scattering of ambient photon 
fields and is also not dependent on any other property of the 
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GC, or that the respective properties are similar to Terzan 5 for 
the six GCs. In this case the X-ray luminosity upper limits only 
have to be compared to the luminosity of Terzan 5, which is 
(2.3±0.3)xlO^^ erg s"^ (see |Eger et al.||2010| assuming again a 
distance of 5.9 kpc). This scenario can be ruled out for every 
GC but NGC 6139, where the luminosity upper limit (UL15) is 
below the expectation, but within the 1 cr statistical uncertainties. 
For this model we also calculated a predicted profile of the radial 
surface brightness (as described in Sect. |3.4| ), which we show as 
stars (red) in Fig.[T] 

4. Conclusions 

Extended and diffuse X-ray emission with a very hard energy 
spectrum has recently been detected from Terzan 5. The aim of 
this work was to search for extended and diffuse X-ray emission 
from other GeV-detected GCs to test whether Terzan 5 stands 
out as an exception or if these properties can be found in a pop- 
ulation of Galactic GCs. We analyzed Chandra data of the six 
GCs M62, NGC 6388, NGC 6541, M28, M80 and NGC 6139 
and did not find any indications for diffuse excess emission on 
similar relative scales as seen from Terzan 5. 

We used the resulting upper limits for the total diffuse X- 
ray luminosity to constrain three different scenarios, based on 
the assumption that the flux measured from Terzan 5 arises from 
the respective emission scenario. In these scenarios diffuse non- 
thermal X-ray emission is either arising from synchrotron ra- 
diation from the relativistic lepton population produced by the 
population of msPSRs, IC emission from already cooled leptons 
accelerated along with the primary hadrons in a GRB remnant, 
or it was simply assumed that every GC hosts a source of diffuse 
X-ray emission with the same luminosity as the one detected in 
Terzan 5. 

We were not able to rule out synchrotron emission in a msP- 
SRs scenario for any of the six GCs. However, this process ap- 
pears to be unlikely for NGC 6388 and M 62. On the other hand, 
a scenario based on a GRB remnant with the same properties as 
the putative GRB remnant in Terzan 5 is ruled out for every GC. 
Lastly, the assumption that each GC hosts a diffuse X-ray source 
with the same luminosity as seen from Terzan 5 is ruled out for 
all GCs but NGC 6139. 

Currently, Terzan 5 is the only GC where a TeV y-ray source 
is found in the vicinity of the cluster. The upper limits on diffuse 
X-ray emission from the other GCs will also be important for 
constraining the emission scenario for TeV y-rays if additional 
GCs are detected in this energy band in the future. 
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